Improved Wilson QCD simulations with light quark masses 
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We present results from simulations using 2 flavours of 0(a)-improved Wilson quarks whose masses 
are about 1/3 of the physical strange quark mass. We present new data on the mass of the singlet 
pseudoscalar meson and evidence of the onset of chiral logarithms in the pion decay constant. The 
previously observed suppression of the topological susceptibility at lighter quark masses is confirmed. 
We report on the performance of the hybrid Monte Carlo algorithm at light quark masses. 
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I. INTRODUCTION 

The first generation of QCD simulations using Wilson 
fcrmions have provided useful information on the hadron 
spectrum but have been restricted to relatively heavy 
quark masses (see for example the annual reviews by 
Aoki and Kaneko Q). Using an improved staggered 
discretisation, the MILC collaboration [3| has probed the 
spectrum with significantly lighter quark masses and very 
promising results are now being obtained 4] from simu- 
lations with light quark masses down to m s /8, a physi- 
cal strange quark mass m s and lattice spacings down to 
0.09 fm. There are theoretical and numerical complica- 
tions associated with this action and the representation of 
lattice flavour symmetries and so simulations with Wilson 
fermions remain an important tool in studying full QCD. 
In due course, simulations with chiral fermions should 
become feasible and provide a further cross-check on the 
results now being obtained. 

In an earlier paper pj the UKQCD collaboration has 
presented lattice QCD results based on two-flavour sim- 
ulations conducted at fixed lattice spacing (a « 0.1 fm) 
and fixed volume around (1.6 fm) 3 . This work also con- 
tained results obtained at fixed gauge coupling {(3 = 5.2) 
also using the standard Wilson gauge action and 0(a)- 
improved Wilson fermions. The lightest sea quark mass 
achieved in these simulations was around 0.6 of the 
strange quark mass (m^/irip w 0.58). There is evidence 
that some hadronic observables suffer finite size effects al- 
ready on these lattice volumes [|| Q • Working at lighter 
quark masses at fixed gauge coupling is expected to en- 
hance further these effects and to provide an even greater 
challenge to the standard Monte Carlo simulation algo- 
rithms. In this paper, we report on the results of an 
attempt to push our analysis towards this lighter quark 
regime in order to (a) uncover any more obvious effects 
of dynamical fermions not hitherto seen with this ac- 
tion U [|| and (b) determine the limits if any of the 
simulation algorithm in its simplest form. A previous 
study had indicated a potential instability at light 
quark masses when using step sizes that were too large 
in the molecular dynamics trajectory part of the update. 



This paper should be seen as an addendum to our pre- 
vious paper in which full details are given of the action, 
simulation methods and basic hadronic measurements. 
Some preliminary results from this work were presented 
in [![l(3. 

Recent related work using improved Wilson fermions 
includes that of the JLQCD collaboration Q\ which uses 
the same action and covers a similar range of quark 
masses to our earlier work but which also includes sim- 
ulations in larger volumes (w (1.8 fm) 3 ). The QCDSF 
collaboration has been conducting simulations comple- 
mentary to those of UKQCD, most recently using an im- 
proved algorithm . Neither of these simulations sets 
has penetrated the lighter quark regime. In each case 
m^/rrip > 0.6. 

The qq+q collaboration has succeeded in simulat- 
ing at one half (or less) of the strange quark mass but 
with rather coarse lattice spacing (a w 0.28 fm). 

The rest of the paper is as follows. In section^ we de- 
scribe the parameters and performance of the new simu- 
lations. The lattice spacing is determined from the static 
potential and decorrelation effects are studied using mea- 
surements of the topological susceptibility. In scction lTTTl 
we present the additions to our previous collection of data 
for the hadron spectrum and meson decay constants and 
use them to search for first signs of chiral logs. Results 
from new measurements of disconnected loops, including 
the r\ mass, are presented in section llVl We present re- 
sults on the topological susceptibility in section Con- 
clusions are drawn in section EU 



II. 



SIMULATIONS WITH IMPROVED WILSON 
FERMIONS 



The simulations were conducted using the standard 
Hybrid Monte Carlo (HMC) algorithm as described 
in [f| using lattice action parameters ((3,c sw ,k) — 
(5.2,2.0171,0.1358) and a lattice volume V = L 3 T = 
16 3 32. A total of 2440 trajectories were accumulated 
at the rate of about 3 per hour on a machine sustain- 
ing 30 Gflops. In physical units, the lattice volume was 
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FIG. 1: Histogram of In AH (positive values of AH only) 
for this ensemble (darker bars) and for a similar sized sample 
from simulations at k — 0.1350. 

L 3 w (1.5 fm) 3 and w 420 MeV before chiral or con- 
tinuum extrapolation. The lattice spacing was estimated 
from the measured hadronic scale parameter |l3j : 

r /a = 5.32(5) . (1) 

The integrated autocorrelation time r of the mean pla- 
quette was found to be 6.9(14) trajectories. This follows 
the previously observed trend of a decrease of r with 
decreasing quark mass. For comparison, t = 16(3) at 
K = 0.1350 where /m p w 0.7. For the present ensem- 
ble, nips/my — 0.44(2) (see later). Although not ex- 
pected, this trend can be accommodated in simple mod- 
els [f| . The mean plaquette was 

(P) = 0.53767(3) . (2) 

These HMC runs are expected to be susceptible to in- 
stabilities (occasional large AH values and zero accep- 
tance) when the fermion force term gets too large (8J. 
We have observed this effect directly in the present sim- 
ulations where we found it necessary to use a step size 
of 1/400. When the step size is this, or smaller, we 
found that 64 bit arithmetic for field storage and matrix- 
vector manipulations was required so as to avoid a se- 
rious loss of acceptance due to rounding errors (see be- 
low). This was true even though we always used full 64- 
bit arithmetic and careful summing techniques for the 
global summation in our AH (energy difference) calcu- 
lations. Fig.[I] shows a histogram of In AH for all trajec- 
tories where AH > 0. For comparison we also show the 
corresponding histogram for the well-behaved simulation 
at (f3, k) = (5.2, 0.1350) described in p|- The sample size 
in each case is the same. 

The anomalous trajectories leading to very large values 
of AH are clearly visible in the simulation at k = 0.1358. 
The simulation was an experimental one and incorpo- 
rated several changes of simulation parameters (step size 
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FIG. 2: The time history of the mean plaquette together with 
the corresponding HMC acceptance (locally averaged over 10 
trajectories). The horizontal bars indicate the precision used 
as discussed in the text. 

and solver accuracy etc.) leading to large changes in ac- 
ceptance. Thus one should bear in mind the possible 
consequences of this on the discussion of autocorrelation 
times and error analysis. 

In Fig.[3we show a time history of the mean plaquette 
along with the average acceptance (integrated over 10 
trajectories). The dramatic drop in acceptance associ- 
ated with a change to 32-bit arithmetic is clearly visible. 
The location of the changes are indicated by the hori- 
zontal bars. The loss of acceptance was not significantly 
dependent on the size of the solver residual in the molecu- 
lar dynamics steps (using a running rather than absolute 
residual) . As noted above, it was primarily dependent on 
the arithmetic used in the matrix-vector calculations. 

Despite the rather chequered history of the configu- 
rations we decided to subject them to physics analysis. 
This ensemble represented an expensive investment in 
computer time and promised to give access to relatively 
light quark masses (by Wilson lattice quark standards). 
In view of the above remarks, one should bear in mind 
the deficiencies of the Markov process that led to their 
generation. We may attempt, in due course, to repeat the 
simulations at this lattice spacing and quark mass using 
improved algorithms, faster machines and larger lattices. 

III. HADRON SPECTRUM AND DECAY 
CONSTANTS 

In this section we report on the light spectroscopy from 
this ensemble. We use a similar analysis to our original 
work on the spectroscopy of nonperturbatively-improved 
clover fermions at (3 — 5.2 Here we concentrate on 
the unitary sector of the theory — with valence quarks 
equal in mass to the sea quarks. We have previously re- 
ported some evidence for chiral logs in a partially 
quenched analysis with K sca = 0.1355 and K V aicncc = 
0.1358. We use fuzzed and local sources and sinks com- 
bined to make a variational fit. To determine the pion 



3 



0.4 



TT 



Meff, ps /c = 0.1358 



0.3 — 



0.2 — 



0.1 




K 


Hadron t- 


-ranj 




ami 


0.1358 


0"+ 


4 


- 15 


0.199(5) 


0.75(11) 


0.1358 


1"- 


4 


- 12 


0.450(14) 1.16(6) 


0.1355 


0"+ 


4 


- 15 


0.282(4) 


0.90(11) 


0.1355 


1— 


4 


- 12 


0.491(7) 


1.17(3) 


0.1350 


0"+ 


4 


- 15 


0.408(2) 


1.19(8) 


0.1350 


1"" 


4 


- 12 


0.585(4) 


1.33(3) 



TABLE I: Masses in lattice units from this calculation. 



a/A 



0.1358 0.0829(26) 0.1457(78) 
0.1355 0.1055(14) 0.1835(44) 
0.1350 0.1336(11) 0.2468(33) 



TABLE II: The raw lattice value of afir is given by using the 
order a improved expression (1 + bAm){afA + CAafp) and we 
tabulate these two contributions. 
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FIG. 3: Effective mass plot for the pseudoscalar channel at 
K = 0.1358. 



decay constant, we fit with two states an order 4 varia- 
tional matrix with fuzzed and local sources using the 75 
and 7475 operators to create a pion. 

To increase statistics, we use quark propagators with 
sources on the time planes 7, 15, and 23 in addition to the 
t = plane for k values 0.1358, 0.1355, and 0.1350. In our 
first published work we only used the quark propagators 
from time plane t = 0. Thus we report new numbers at k 
= 0.1350 and 0.1355 for 7r and p with reduced error bars. 

In Fig.|3|we show the effective mass plot for the pseu- 
doscalar channel (PS) at k — 0.1358. The results from 
the fits are in Tables [I] [HI We investigated the stability 
of the fits in a number of ways. Since the spectrum fits 
used two states, the ground state is expected (and in- 
deed found) to be very stable while the excited state val- 
ues should be regarded as indicative, although these val- 
ues are consistent with the expected lightest multi-body 
states at 3m,r and 2rn 7r (fc = 2ir/L), for pseudoscalar and 
vector, respectively. 

Our values of m^L « 4, so we should expect some finite 
volume effects. Using a similar formalism and parame- 
ters (the only difference being to use csw = 2.02 rather 
than 2.017), JLQCD |(j have explored this for the hadron 
spectrum using L — 12, 16, 20. Here we discuss this at 
their lightest K-value (0.1355), where they do see evidence 
of a finite size effect (masses lower at larger volume) al- 
though it is not very significant statistically for the two 
larger volumes. We have a larger number of trajectories 
at L = 16 than JLQCD, and there are some statistically 
significant differences between our results |5j] and theirs 



for that volume. However, for the pseudoscalar and vec- 
tor mesons, the differences between our new determina- 
tions at L = 16 (given above) and theirs at K = 0.1355 
are not statistically significant. Using our newer results 
to make the comparison with L = 20 and L = 16, then 
suggests that there is a decrease of the pion mass as the 
volume increases by about 4 ± 2%. 

We have also computed the pion decay constant for 
this data set. It has been clear from the work of the 
GF11 group that the unquenching effects are larger in 
decay constants than for masses 15j, although the sys- 
tematic errors on decay constants can be large due to 
truncation of perturbative series. The ratio fx/ fn (for 
which renormalisation factors should largely cancel) is 
underestimated in quenched lattice QCD [l(j ■ 

A critical goal of lattice gauge theory calculations is to 
detect the presence of chiral logs in observables. The loop 
corrections from chiral perturbation are non-analytic in 
the parameters of the lowest order chiral Lagrangian, and 
hence provide a good check that the lattice calculation is 
in the regime where chiral perturbation theory is valid. 
A particularly appropriate way to look for these chiral 
logarithms is in the pseudoscalar decay constant. 

We are working, because of computational constraints, 
at finite lattice spacing. The formalism of chiral pertur- 
bation theory can be extended to cover this case [13, [l8| , 
but at the cost of additional parameters. We choose in- 
stead to make a comparison with the continuum predic- 
tions of chiral perturbation theory. 

Our results for are shown in Fig.0] We have ex- 
tracted the values from fits as discussed above and then 
applied the rotations and corrections appropriate for an 
order a improved formalism. We have used the same per- 
turbative formulation of these corrections (and also the 
same prescription for Za ) as employed by JLQCD, in 
order to facilitate comparisons. Moreover, since they use 
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FIG. 4: The pseudo-scalar decay constant in units of ro from 
UKQCD and JLQCD versus k. 



a different prescription for evaluating r , we have applied 
our determination of ro to their data. The comparison 
is shown in Fig.^J The agreement at our two common k 
values is adequate and we present here a new determina- 
tion at a lighter k value. The significant feature of this 
new result is that it shows a curvature versus m q . 

This curvature may be the first evidence in this study 
of the chiral logarithm at work. To explore this we com- 
pare our result with some continuum chiral models for 
which the chiral logarithm has a fixed coefficient (given 
by /„.). For Nj — 2 flavours of degenerate quark, the 
standard chiral perturbation theory result to one loop is: 
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FIG. 5: The pseudoscalar decay constant in units of ro from 
UKQCD versus the squared pseudoscalar meson mass. Also 
shown is an expression including chiral perturbation theory 
terms to order m 4 which has been fitted (see ref. |19| where 
we use /i = 0.75 GeV and ri?(/i) = —2) to agree with the 
experimental values of f n and /k which are shown (*). An 
estimate [Tsll of the finite size effect expected from chiral per- 
turbation theory (to order m 2 ) is shown by the vertical lines. 
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This expression has an unsatisfactory behaviour at 
large m, where chiral perturbation theory should not ap- 
ply anyway. Moreover even at the K mass, the quar- 
tic terms in a chiral perturbative treatment are signif- 
icant [lj|. We can thus either concentrate on the cur- 
vature implied at small m or modify the expression phc- 
nomenologically. We illustrate this behaviour in Fig.[5]by 
using an empirical determination |l9j | of the terms arising 
in chiral perturbation theory up to m . This shows the 
curvature to be expected at small pseudoscalar mass in 
a large volume. Unfortunately the overlap between our 
data and the region of validity of the chiral approach to 
this order is small. 

As discussed above, finite volume effects should be 
important, since for this lightest quark mass, we have 
m^L — 3.2. As has long been known, continuum chiral 
perturbation theory in a finite volume gives ex plic it pre- 
dictions. This has been explored theoretically [l!| using 
a range of different treatments of the chiral formalism, 
showing that we are in the region where chiral models 
do give rise to significant effects. Moreover the finite size 
effects arise from chiral loops and so are from the same 



source as the logarithmic corrections. Using L = 1.5 fm 
and ttl-jy — 400 MeV, which are close to our values, a 
range of different levels of approximation yields a 
relative increase of the finite-volume pion mass over the 
infinite volume mass of 0.6 to 2%. This can be compared 
to the shift of 4 ± 2% obtained from lattice comparisons, 
as discussed above. This is not a statistically significant 
comparison, but it does indicate that the magnitude of 
finite size effects to be expected theoretically is consistent 
with that seen on the lattice. 

In the Gasser-Leutwyler chiral approach, the relative 
finite size effects for / w will be four times as large as 
those for the pion mass, and of opposite sign. Again, 
this finite size effect comes from the one loop term that 
gives the chiral logarithm. We show the effect of this 
shift for the Gasser-Leutwyler approach 0, 12(1 l2ll [2^ 
(corresponding to 3.4% for the above values of L = 1.5 fm 
and m = 400 MeV) in Fig.|SJ Note that the curvature we 
observe is equivalent to a decrease of f n at k — 0.1358 of 
about 8%, which could thus be ascribed entirely to finite 
size effects only if the one loop estimate was less than 50 
% of the total. Since the finite size estimates come from 
chiral models, they will also include a chiral logarithm 
which will yield curvature. 
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We conclude that the finite volume effects will enhance 
the curvature at small pion masses, as indeed we find. 
Hence our results are in qualitative agreement with chiral 
perturbation theory. 

We note that our data suggests an extrapolation to the 
chiral limit which would give a value for f n below the ex- 
perimental value of 131 MeV. Since the perturbative cor- 
rection to Za is of order 25% at first order, we expect pos- 
sible systematic errors of up to 5% from the next order, 
which is only estimated by using tadpole- improved meth- 
ods. This error budget for Za is confirmed by results 
from quenched studies where the non-perturbative eval- 
uation of Za gave a value 4% different from the tadpole- 
improved one loop result used here. We chose to use a 
value of rp = 0.525 fm to set the scale following previous 
work 0,123 an d this value is uncertain to 5%. There are 
further errors in the lattice determinations coming from 
having a sea quark mass which is too large, from neglect 
of the strange sea and from finite a effects (as well as 
the finite volume effects we have discussed above). This 
covers the apparent discrepancy of 10% seen in the fig- 
ure comparing our result to the experimental values of 
/„■ and f K . 

Although it may appear that reproducing f n and fx 
from lattice QCD has no immediate experimental impact, 
that is not quite true. The search for chiral logs in decay 
constants is currently one of the most important topics 
in heavy- light physics [24l l2fij . The error on the ratio 
of the fs 3 1 fs has recently been increased, because the 
chiral log term has not been observed in lattice data (2|| . 
The ratio /s a //s is an important QCD quantity for the 
unitarity checks of the CKM matrix. It will become more 
important once B s mixing is measured at run II of CDF. 
As has been noted by many authors [27], , the chiral 
log structure of / ff and fg are rather similar. Hence, a 
detection of chiral logs in f v is an indication that the pa- 
rameters of the unquenched calculation are close to where 
chiral logs may occur in the heavy-light decay constant. 

Evidence for chiral logs in both heavy and light QCD 
has been claimed in unquenched calculations with im- 
proved staggered quarks 29] . Particularly because of the 
complexity of the chiral perturbation theory calculations 
for staggered fermions j33, US] j we feel a cross check on 
the improved staggered calculations is very valuable, even 
if ultimately unquenched clover fermions do not allow us 
to control all the systematic errors such as lattice spacing 
dependence. 



IV. THE MASS OF THE SINGLET 
PSEUDOSCALAR 

The large splitting between the mass of the rj and the 
octet of light pseudoscalars is thought to be caused by the 
complex vacuum of QCD and the anomaly in the axial 
current. There is a lot of activity in trying to reconcile 
the mechanism behind the mass splitting between the 
mass of the rj and the masses of the octet. In particular 
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FIG. 6: The ratio of disconnected D to connected C contribu- 
tions to the SP two-point correlator (C — D) at k = 0.1358. 
The operators used for the pion are local (L) of fuzzed (F) 
and either 75 (P) or 7574 (A). 



the questions raised by Witten about the consistency of 
the solution of the U(l) via instantons with the large N c 
limit [32j is topical. There are also many phenomenolog- 
ical puzzles in which the rj is involved. See the review 
by Bass |33j for a review of experiments with rj and rj 
as decay products. A first principles calculation of the 
structure of the rj would be helpful. 

In the real world the mass of the rj is also determined 
by the mixing between the singlet and octet mesons. 
We introduce the notation: NP is the nonsinglet pseu- 
doscalar and SP is the singlet pseudoscalar. This mix- 
ing can be estimated from partially quenched two flavour 
QCD 0. 

We use a similar methodology to that used in a ear- 
lier study on a coarser lattice |34|. The fermion loop 
was computed using complex Z2 noise using the "two- 
source" trick we used in the calculation of the non-singlet 
scalar 35]. We used 100 noise samples. Using fuzzed 
smearing functions as a basis, we fitted to a 2 by 2 matrix 
of correlators using "factorising fits" . The A4 operator 
also couples to the pion so we sometimes use a basis of 
4 smearing functions. To show the quality of data we 
plot the ratio of disconnected to connected correlator in 
Fig .0 Note that unitarity requires that the SP correla- 
tor (C — D) is positive so that D/C < 1, as we do indeed 
find. 

The results from the fits are in Tabl eHTTl As expected, 
the fits with a larger basis of smearing functions have 
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TABLE III: Fits to the SP particle. Correlators LL, FL and 
FF are used in each case with a) having PP only; b) having 
PP, AP and AA, while c) has all of PP, AP, PA and AA, 
where P is the pseudoscalar coupling (75) and A is the time 
component of the axial (7574). The symbol 1 implies momen- 
tum^. 



smaller error bars. 

For the fits to the k = 0.1355 data we could obtain a fit 
with three exponentials. We regard the last exponential 
as representing the truncation error, hence we have in- 
formation on one excited state. The mass of the excited 
state at k — 0.1355 is 1.5(2) GeV with unknown sys- 
tematic errors from the lack of continuum extrapolation 
and chiral extrapolation. This is encouragingly close to 
the the mass of the ?y(1295) and 77(1440) mesons. When 
the systematics in the lattice calculation are under con- 
trol, the comparison with experiment will also require an 
understanding of the mixing. 

For the k = 0.1358 data, we are unable to obtain con- 
vincing and consistent fits. Those shown illustrate the 
problem. We use a different method to show the impact 
of the data: namely a direct comparison between the SP 
correlators at the two kappa values, see Fig.[7| For the LL 
and FF correlators (with 75 at source and sink) the ra- 
tios are consistent with amsp(0. 1358) — amsp(0. 1355) = 
—0.03(3) as shown, where we have relied more on the FF 
data since it has a larger contribution from the ground 
state. This mass difference can then be used to ex- 
tract an estimate of the SP mass at K = 0.1358 of 
amsp — 0.40(5). Note that although the SP mass is ap- 
proximately constant as the quark mass is reduced, the 
pion mass (in lattice units) decreases by about 0.09 and 
thus the difference between the SP mass and the pion 
mass increases. This large mass splitting is consistent 
with the steep rise shown in Fig-El 

Because the signal to noise is so poor, we also explore 
singlet correlators with non-zero momentum. One has to 
be careful since, at non-zero momentum, the axial (A4) 
pion operator (35(74) has contamination from a\ — so fac- 
torising fits need additional care for the A4A4 term. For 
K = 0.1355, we can fit for momentum (1,0,0) (in units 
of 2-k/L). For k = 0.1358, we find no useful additional 
constraint from the momentum non-zero correlators. 

The mass of the SP meson in two flavour QCD is 
not immediately available from experiment since the well 
known mixing between the rj and rj obscures this is- 
sue. By assuming some mixing scheme, we can obtain 
an estimate of the mass of the SP particle. In our 




FIG. 7: The ratio of SP (Nf = 2 r\ meson) two point correla- 
tions from k = 0.1358 to k — 0.1355. The curve is described 
in the text. 



previous analysis 



, the mass of the SP meson was 
Using values consistent with phe- 



m o = V 771 ss + % x 
nomenology and our previous lattice data (m ss = 0.695 
GeV, x ss = 0.13 GeV 2 ), we obtained 0.861 GeV as the 
mass of the SP meson in Nf = 2 QCD. 

Another approach is to use the Witten-Veneziano ex- 
pression. The SESAM/TxL collaboration ^36| obtain 715 
MeV for mass of the SP in N f = 2 QCD. 

The chiral extrapolation formulae used in the lattice 
QCD literature j^E have either used the mass or 
square of the mass of the SP meson linear in the quark 
mass. It would be clearly better to have a more theoret- 
ical justification of the light quark mass dependence, al- 
though "traditional" chiral perturbation theory is not re- 
liable at mass scales appropriate to the SP state [3Sl l39]| . 

In Fig.|Hl the world data for the SP mass is plotted 
as a function of the pseudoscalar to vector mass ratio. 
Our point at k = 0.1355 is consistent with the data from 
other groups. The mass at k = 0.1358 lies above the 
trend from larger quark mass, but we do expect the SP 
mass to go to a non-zero constant as the quark mass van- 
ishes. Indeed our semi-phenomenological estimate given 
above was that the SP mass in the chiral limit is 861 
MeV. This value is very consistent with the flattening 
behaviour indicated by our data point. 
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FIG. 8: World eta data from lattice studies with two degen- 
erate sea quarks. The bursts and squares are the values from 
SESAM from the truncated eigenvalue analysis |4f| and Z2 
noise measurements [3(J respectively. The crosses are from 
UKQCD (k = 0.1355 and k = 0.1358 from this work and 
K = 0.1398 from ref. 0)- Results from CP-PACS ^ are 
shown by a diamond at /3 = 2.1, an octagon at j3 = 1.95 and 
a fancy plus at j3 — 1.8 



For their final result, CP-PACS 37] quoted the mass of 
the SP particle as 0.960(87) tg;!^ GeV. This high value 
arises from the continuum extrapolation. The central 
value for the mass was obtained by linearly extrapolated 
in lattice spacing with a x 2 /dof of 4.2. A quadratic 
extrapolation in the lattice spacing had a x 2 /doi of 
2.8, with the resulting mass of 0.819(50) GeV. The use 
of a linear extrapolation with lattice spacing is consis- 
tent with the rest of the spectroscopy program of CP- 
PACS [4f| . The large errors in the mass of the SP meson 
from CP-PACS represent the variation from the different 
continuum extrapolations. 



V. TOPOLOGY 

The topological charge, Q, and its associated suscepti- 
bility, 



X 



L 3 T 



(4) 



are expected to be especially sensitive to the presence and 
properties of the sea quarks in QCD. As the mass of the 
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FIG. 9: A time history of the topological charge, Q, with 
a histogram using unit-sized bins, and the Gaussian curve 
derived from the central value of topological susceptibility 
(and one standard deviation either side as outlying curves). 



sea quarks is reduced (towards the chiral limit), the topo- 
logical susceptibility is suppressed below its quenched 
value. For sufficiently lar ge v olumes, the leading order 
chiral behaviour is |42J, |43j, |44j 



X(ml) = 



(/•7T«V) 5 



(5) 



for two degenerate flavours (and using a normalisation 
where f n ~ 93 MeV). As m q — » oo, x —* X qn the 
quenched value which is around (180 MeV) 4 . The higher 
order corrections in eq. JSJ must therefore introduce a 
neg ative curvature at some intermediate quark mass 

On the lattice, the topological susceptibility becomes 
renormalised relative to the continuum value, Xcont, both 
multiplicatively and additively |47]| : 



X 



Z 2 Xcont 



M 



(6) 



Broadly speaking, M > arises from the presence of 
"dislocations" : short range fluctuations in the gauge field 
that masquerade as small instantons. Z < 1 reflects the 
breaking of scale invariance on the lattice, whereby small 
instantons have a topological charge less than unity. At 
large quark masses (or in the quenched theory) the first 
term in eq. © dominates, suppressing \ a t non-zero a. 
In the chiral limit, however, M dominates and is non-zero 
even after smoothing. The topological susceptibility then 
shows strong discretisation effects that act to increase \ 
at finite lattice spacing. 

Whilst comparing "matched" ensembles at fixed lat- 
tice spacing will control discretisation effects, these oppo- 
site trends imply that they will not cancel away entirely. 
The net effect is that is that any chiral suppression of 
the topological susceptibility at given a relative to the 
quenched value at an equivalent lattice spacing will be 
less than in the continuum limit. 

We measure Q using the method of 45J: n c = 10 cool- 
ing sweeps are applied using the Wilson gauge action. 



8 



0.08 



0.06- 




FIG. 10: The topological susceptibility as a lunction ol the 
lightest pseudoscalar ("pion", « V alence = K sca ) mass for two 
flavours of C(a)-improved fermions. Quenched values for this 
range of lattice spacing are shown as a bar on the RHS of the 
plot. 



[in (rom-rr) 2 ] fit, or an interpolating fit across the chi- 
ral range [5(j], we see that the slope near the origin of 
Fig. ^3 is slightly greater than that expected from the 
continuum value of f^. As discretisation effects are ex- 
pected to increase the pseudoscalar decay constant, this 
is in agreement with theoretical expectations. In a forth- 
coming paper discretisation effects will be examined in 
more detail |52| . 

Finally, it is interesting to compare these results with 
those obtained for three flavours of improved staggered 
sea quarks ^7|- For similar lattice spacings, a ~ 0.09 fm, 
the topological susceptibility was non-zero and roughly 
constant below (rom^) 2 ~ 2, when x presumably became 
dominated by M. No statistically significant evidence 
for such a trend can yet be seen for the C(a)-improved 
action, however. 



VI. CONCLUSIONS 



Ten cools strikes a good balance between adequate sup- 
pression of these ultraviolet dislocations and excessive de- 
struction of the long range topological structure [45j. A 
reflection-symmetrised "twisted plaquette" lattice topo- 
logical charge operator is used [43 ■ 

As discussed previously, good decorrelation of Q is seen 
in the simulation, and the histogram of the populations 
of the different topological sectors in Fig.[§] has the ex- 
pected Gaussian form. We find (Q) = —0.33(29), consis- 
tent with zero. The susceptibility is x = 0.292(45) x 10~ 4 , 
or 0.284(34) x 10~ 4 if we subtract terms in (Q). We plot 
the latter result as the leftmost point in Fig.MiJJ along- 
side previously published UKQCD results [M l45 l liM lia 
IHcl IFlf . We show also the equivalent quenched suscepti- 
bility as a bar whose length reflects both the statistical 
uncertainty in the quenched measurements and the small 
variation in lattice spacing across the ensembles depicted, 
despite the presence of discretisation effects. We see very 
clear evidence for a strong chiral suppression of x relative 
to the quenched theory, driven by the sea quarks. 

For eq. JSJ) to hold, we require xls = (/ttTOtt) 2 ^ ^> 1 
|44|. Using the continuum value of f n (which is lower 
than the value at finite lattice spacing) we find xls =11- 
Finite volume effects appear to be less significant for x 
than, for instance, the light hadron spectrum. Studies for 
this action indicate that ils ^ 10 is sufficient for such 
finite size effects to be within the statistical uncertainty 
in X |H. 

As discussed above, comparison of data at finite lat- 
tice spacings with continuum predictions must be made 
cautiously. That being said, Fig. is very encourag- 
ing and is evidence for the improved chiral properties of 
the C(a)-improved action. Performing a leading order 



Although we are excited to have finally reached a re- 
gion of parameter space where unqucnched clover calcu- 
lations are starting to see chiral logs and uncover inter- 
esting behaviour in the mass of the singlet pseudoscalar 
meson, it is not yet clear how to improve on these results. 

If we im ple mented some of the new updating algo- 
rithms |53l l54| for clover fermions, in principle we could 
work at a larger volume with fewer of the problems re- 
ported in section [H] However, concerns have been raised 
about the interaction of the Wilson gauge action 

with the clover fermion action in unquenched calcula- 
tions. Once improved gauge actions have been incor- 
porated into the non-perturbative clover improvement 
formalism for fermions, we could study the chiral log 
structure, at fixed lattice spacing 0, This would 
allow important comparisons with the results from other 
unquenched lattice QCD calculations that use different 
fermion formalisms 41 . 



Acknowledgments 

We acknowledge the support of the U.K. Parti- 
cle Physics and Astronomy Research Council un- 
der grants GR/L22744, GR/L29927, GR/L56374, 
PPA/G/O/1998/00621, PPA/G/O/1998/00518, 
PPA/G/S/1999/00532, PPA/G/O/2000/00465 and 
PPA/G/O/2001/00019. A.H. is supported by the U.K. 
Royal Society. We thank EPCC for time on Lomond. 
We are grateful to the ULgrid project of the University 
of Liverpool for computer time. The authors acknowl- 
edge support from EU grant HPRN-CT-2000-00145 
Hadrons/LatticeQCD . 



[1] S. Aoki, Nuc l. Phys. Proc. Suppl. 94, 3 (2001), [2] T. Kaneko, Nucl. Phys. Proc. Suppl. 106, 133 (2002), 
|hep-lat/00li074l 



9 



|hep-lat/01lT005| [28 
[3] C. W. Bernard e t al, Phys. Rev. D64, 054506 (2001), 

hep-lat/0104002 [29 
[4] HPQCD, C. T H. Davies et al, Phys. Rev. Lett. 92, 

022001 (2004), hep-lat/0304004 [30 
[5] UKQC D, C. R AUton e t al, Phys. Rev. D65, 054502 

(2002) , hep-lat/0107021 [31 
[6] JLQCD, S. Aoki et ar7"Phys. Rev. D68, 054502 (2003), 

hep-lat/0212039 [32 
[7] QCDSF, A. Ali Khan et al, Nu cl. Phys. Proc. Suppl. [33 

119, 419 (2003), hep-lat/0209111 
[8] UKQCD, B. Joo et al, Phys. Rev. D62, 114501 (2000), [34 

hep-lat/0005023 

[9] UKQCD, ATTT Irving, Nuc l. Phys. Proc. Suppl. 119, [35 

341 (2003), hep-lat/0208065 
[10] UKQCD, D. J. Hepburn, Nucl. Phys. Proc. Suppl. 119, [36 

338 (2003), hep-lat/0209036 
[11] QCDS F, A' Ali Khan T Tal, Phys. Lett. B564, 235 [37 

(2003) , [hep4a t /0303026 
[12] qq+q, F. Farchioni, C. Gebert, I. Montvay, E. Scholz, [38 

and L. Scorzato , Phys. Lett. B561, 102 (2003), 

hep-lat/0302011 [39 
[13] R. Sommer, Nucl. Phys. B411, 839 (1994), 

hep-lat/9310022 [40 
[14] UKQCD, A. C. Irving, C. McNeile, C. Michael, K. J. 

Sharkey, and H. Wittig, Phys. Lett. B518. 243 (2001). [41 

hep-lat/0107023 

[15] F. Butler, H. Chen, J. Sexton, A. Vaccarino, and [42 
D. Weingarten, Nucl. Phys. B421, 217 (1994), [43 
hep-lat/9310009 

[16] ALPHA, J. Heitger, R. Sommer, and H. Wittig, Nucl. [44 

Phys. B588, 377 (2000), hep-lat/0006026 
[17] S. Aoki, P hys. Rev. D68, 054508 (2003), [45 

hep-lat/0306027 

[18] O. Bar, G. R upak, and N. Shoresh, (2003), [46 

hep-lat/0306021 [47 

[19] G. Colangelo and S. Durr, (2003), hep-lat/0311023 
[20] J. Bijnens, G. Colangelo, G. Ecker, J. Gasser, and M. E. [ 

Sainio, Phys. Lett. B374, 210 (1996), hep-ph/9511397 
[21] J. Bijnens, G. Colangelo, G. Ecker, J. Gasser, and M. E. [49 

Sainio, Nucl. Phys. B508, 263 (1997), hep-ph/9707291 
[22] J. Bijnens, G. Colangelo, and P. Talavera, JHEP 05, 014 [50 

(1998), hep-ph/9805389 
[23] UKQCD, A. Dougall, R. D. Kenway, C. M. May- [51 

nard, and C. McNeile, Phys. Lett. B569, 41 (2003). 

hep-lat/0307001 [52 
[24] N. Yamada, Nu cl. Phys. Proc. Suppl. 119, 93 (2003), [53 

hep-lat/0210035 

[25] A. S. Kronfeld, (2003), hep-lat/0310063 [54 
[26] A. S. Kronfeld and S. M. Ryan, Phys. Lett. B543, 59 [55 

(2002), hep-ph/0206058 [56 
[27] MILC, O Bernard et al, Phys. Rev. D66, 094501 (2002), 

hep-lat/0206016 



Phys. Rev. D67, 074503 
Eur. Phys. J. C17, 623 

Phys. Rev. D63, 074019 
Phys. Rev. D65, 054505 



D. Becirevic, S. Fajfer, S. Prelovsek, and J. Zupan, Phys. 
Lett. B563, 150 (2003), |hep^ph/021127l| 

M. Wingate, C. T. H. Davies, A. Gray, G. P. Lepage, and 
J. Shigemitsu, (2003), |hep^ph/0311130| 
C. Aubin and C. Bernard, Phys. Rev. D68, 034014 
(2003), hep-lat/0304014 

C. Aubin and C. Bernard, Phys. Rev. D68, 074011 
(2003), hep-lat/0306026 

E. Witten, Nucl. Phys. B149, 285 (1979). 

S. D. Bass, Phys. Scripta T99, 96 (2002), 
hep-ph/0111180 

UKQCD, C. Mc Neile and C Michael, Phys. Lett. B491, 
123 (2000), hep-lat/0006020 

UKQCD, C. McNeile and TJ. Mi chael, Phys. Rev. D63, 
114503 (2001), hep-lat/0010019 

TxL, T. Struckmann et al, Phys. Rev. D63, 074503 
(2001), hep-lat/0010005 
CP-PA CST^V- L Lesk et al, 
(2003), hep-lat/0211040 
R. Kaiser and H. Leutwyler, 

(2000) , hep-ph/0007101 
B. Borasoy and S. Wetzel, 

(2001) , hep-ph/0105132 
CP-PACS, ~A. Ali Khan et al. 

(2002) , hep-lat/0105015 

H. Neff, T. Lippert, J. W. Negele, and K. Schilling, Nucl. 

Phys. Proc. Suppl. 119, 251 (2003), hep-lat/0209117 

R. J. Crewther, Phys. Lett. B70, 349 (1977). 

P. D. Vecchia and G. Veneziano, Nucl. Phys. B171, 253 

(1980). 

H. Leutwyler and A. Smilga, Phys. Rev. D46, 5607 
(1992). 

UKQCD, A. Hart and M. Teper, (2000), 
hep-ph/0004180 

S. Durr, Nucl. Phys. B611, 281 (2001), |hep4at70 103011 
MILC, C. Bernard et al, Phys. Rev. D68, 114501 (2003), 
hep-lat/0308019 

UKQCD, A. Hart a nd M. Teper, Nuc l. Phys. Proc. Suppl. 

83-84, 476 (2000), hep-lat/9909072 

UKQC D, A. Hart, Nuc l. Phys. Proc. Suppl. 106, 575 

(2002) , hep-lat/01 10167 

UKQC DTX Hart and M . Teper, Phys. Lett. B523, 280 
(2001), hep-lat/0108006 

UKQCD, A. Hart and M. Teper, (2000), 

hep-lat/0009008 

UKQCD-QCDSF, A. Hart, (2004), hep-lat/0401015 
M. Hasenbusch and K. Jansen, Nucl. Phys. B659, 299 

(2003) , |hep-lat/02lT042l^ 

M. Hasenbusch, (2 003),|hep4a t/0310029 

K. Jansen, (2003) , |hep-lat /0311039| 

ALPHA, R. Sommer et al, (2003), hep-lat/0309m] 



